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Abstract 

The pion and kaon electromagnetic form factors Fm{Q'^) are calculated 
at the leading order of pQCD using the running coupling constant method. 
In calculations a dependence of the mesons distribution amplitudes on the 
hard scale is taken into account. The Borel transform and resummed 
expression for Fm{Q'^) are found. 
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1. Investigation of the infrared (ir) renormalon effects in various inclusive and 
exclusive processes is one of the important and interesting problems in the perturba- 
tive QCD (pQCD) [1]. It is known that ir renormalons are responsible for factorial 
growth of coefficients in perturbative series for the physical quantities. But these 
divergent series can be resummed by means of the Borel transformation [2] and 
the principal value prescription [3] and effects of ir renormalons can be taken into 
account by scale-setting procedure asiQ"^) (^si^^Pif {Q'^))Q'^) the one-loop 
order results. Technically, all-order resummation of ir renormalons corresponds to 
the calculation of the one-loop Feynman diagrams with the running coupling con- 
stant as{—k'^) at the vertices or, alternatively, to calculation of the same diagrams 
with non-zero gluon mass. Both these approaches are generalization of the Brodsky, 
Lepage and Mackenzie (BLM) scale-setting method [4] and amount to absorbing cer- 
tain vacuum polarization corrections appearing at higher-order calculations into the 
one-loop QCD coupling constant. Studies of ir renormalon problems have opened 
also new prospects for evaluation of power suppressed (higher twist) corrections to 
processes' characteristics [5]. 

Unlike inclusive processes exclusive ones have additional source of ir renor- 
malon contributions [6-9]. Thus, integration over longitudinal fractional momenta 
of hadron constituents in the expression of the electromagnetic (elm) form factor 
[6-8] or, in the amplitude of exclusive subprocesses [9] generates ir renormalon cor- 
rections. 

In this work we calculate the pion and kaon elm form factors using the running 
coupling constant method and by taking into account dependence of the mesons 
distribution amplitudes on the hard scale Q^. 

2. In the context of pQCD the meson elm form factor has the form, 

FniQ') = C C dxdy<Pl,{y, Q')Th{x, y; Q\ as{Qi'))4>M{x, Q^), (1) 
Jo Jo 

where — —q^ is the square of the virtual photon's four-momentum. In Eq.(l) 
(pAtix, Q^) is the meson distribution amplitude, containing all non-perturbative 
hadronic binding effects, whereas TH{x,y;Q'^,as{Q'^)) is the hard-scattering am- 
plitude of the subprocess qq' + 7* — > qq' and can be found using pQCD. At the 
leading order Th is given by the following expression 

2 as{Q\l-x){l-y)) l asjQ^xy^ 
3 {l-x){l-y) 3 xy J " ^ ^ 

In Eq.(2) is taken as the square of the momentum tranfer of the exchanged 
hard gluon in corresponding Feynman diagrams, Cp = 4/3 is the color factor. 

One of the important ingredients of our study is the choice of the meson distri- 
bution amplitude (pM^x^Q"^). In the framework of pQCD it is possible to predict 
the dependence of 4>m{x, Q^) on using evolution equation, but not its shape (its 



TH{x,y;Q',as{Q')) 
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dependence on x). An information about the shape of the distribution amphtude 
4>m{x, Q^) can be obtained by means of the QCD sum rules method or other non- 
perturbative approaches and as a result some model functions for (f)M{x,Q'^) niay 
be proposed. Model distribution amplitudes for the mesons are, therefore, objects 
of contradictory statements and conclusions made in the literature. Such discus- 
sions is particularly intensive in the case of the pion amplitude (f).,r{x,Q^) [10,11]. 
Nevertheless, in this article for the pion we take "traditional" model distribution 
amplitudes proposed in the literature [12-14]. They have the following form 



where <Plsy{x) is the pion asymptotic distribution amplitude. 



(f^lsyi^) = V3Ux{l-x), 



(3) 



(4) 



and f^^ is the pion decay constant; — 0.093 GeV. 

The constants a, b, c in Eq. (3) were found by means of the QCD sum rules 
method at the normalization point /iq — 0.5 GeV and have the values 



a — 0, b — 5, c — 0, CZ amplitude [12], 
a — —0.1821, b = 5.91, c = 0, amplitude from Ref .[13], 
a = 0.6016, 6 =-4.659, c= 15.52, amplitude from Ref .[U]. 



(5) 



The meson distribution amplitude evolves in accordance with the expression 
[12,15], 



n=0 



as {Q') 

«5 il^o) 



7n//3o 
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In Eq.(6) ^Cf/'^{2x — 1)| are the Gegenbauer polynomials, (3q = 11 — 2n//3 is the 
QCD beta-function's first coefficient and 7„ is the anomalous dimension. 



In 



4 
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n+l 1 

1 = + 4y - 

(n + l)(n + 2) f^,j 



The pion distribution amplitude (Eq.(3)) can be rewritten in the form (6). But 
after defining of the coefficients and taking into account the evolution of 4>n{x, Q"^) 
on Q^, for our purposes it is instructive to give to the distribution amplitude its old 
form, namely 



K{x,Q^) 



^asy 



[x) a + (5{2x-lf + '^{2x-iy- 



(7) 
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In Eq.(7) the new coefficients a, are functions of 

-A2 



a 



a + 



b Uc 
5 ^ 105 



— (1 - ^4) , 
21^ ^' 



14c 



Here An is 



(3^bA2 + — {A,-A,), 



7 = cAi. 

asm 



(8) 
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Using tlie same metliod for tfie kaon wc find 

(f>K{x, g') = (pasvi^) fa + 7 (2x - 1) + /3(2x - 1)' + 5{2x - if 



a 



'asy\ 

a+^(l-^2),/3 = 6^2,7 
o 



(9) 
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(^1 - ^3) , 5 = cA^. 



In Eq.(9) (j)asy{^) same form as in Eq.(4), but with replaced by fx — 0.122 

GeV and constants a,b,c take values [12] 

a = 0.4, 6 = 3, c = 1.25. 

The QCD running coupling constant asiQ"^) in Eq.(2) suffers from ir singular- 
ities associated with the behaviour of the as{Q'^) in the soft regions a; — > 0, y — > 
0; X — > 1, y — > 1. Therefore, Fm{Q'^) can be found after proper regularization of 
CKsiQ^) in these soft end-point regions. For these purposes let us relate the run- 
ning coupling as{XQ^) in terms of cks{Q'^) by means of the renormalization group 
equation 



asiQ' 



1 + (Q;5(Q2)/3o/47r) In A" 



(10) 



where as{Q'^) is the one-loop QCD coupling constant. 

3. As was demonstrated in our previous works [6-8], integration in Eq.(l) using 
Eqs.(2),(10) generates ir divergences and as a result for Fm{Q'^) wc get a pcrturba- 
tive series with factorially growing coefficients. This series can be resummed using 
the Borel transformation [2], 



(167r/. 



M 



duexp 



Anu 



B 



Q^Fm {u), (11) 



where B[Q'^Fm\{u) is the Borel transform of the corresponding perturbative series 
[7,8]. 



4 



But simple trick allow us from Eq.(l) directly obtain Eq.(ll). In fact, after 
changing in Eq.(l) the variables x and y io z = ln(l — x) an w = ln(l — y) (or 
z — lax, w = h\y) and using formula 



f 

Jo 
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e-"(*+"'+^)dii = , t 



t + w + z' PoOisiQ'^ 



and after integration over z, w we find for the meson elm form factor the expression 
(11), where the Borel transform has the form 



N 



B 



m„(g2) i„(g2y 



Q^i^M N = E + ■ (12) 



The exact expressions for m„((5^), In(Q^) are given in the Appendix. 

The Borel transform (12) has double and single poles at u — n. Then the 
resummed expression (11) can be calculated with the help of the principal value 
prescription [3], 



+ (1„ + m„ In A) 



1 \-n 1 n — - V \r, 

n A^ 



(13) 



/^O n=l 

where Li (A) is the logarithmic integral [16], 

U{\) = P.V. A = Q7A^ (14) 

JO mx 

In Eq.(13) we have taken into account the dependence of the distribution am- 
plitude (I)m{x, Q'^) on the scale and it is the generalization of our results for the 
pion and kaon elm form factor. Indeed, if we switch off this dependence {A^ = 1), 
then Eq.(13) coincides for a = 0, 6 = 5, c — and N — A with our result for the 
pion Rcf.[7] and for a = 0.4, b — 3, c — 1.25 and N — 5 with the kaon elm form 
factor obtained in Ref. [8]. 

4. For studing phenomenological consequences of the ir renormalon contribu- 
tions to Q"^ Fm^Q"^) 1 it is instructive to introduce the ratio Rm = [Q'^Fm{Q'^)Y^'^ / 
[(5^Fm(Q^)]°, where [(5^-Fm(Q^)]° is the meson elm form factor calculated in the 
context of the frozen coupling approximation. 

In the frozen coupling approximation, for example, the pion form factor is 



In Fig.l the dependence of the ratio on (CZ distribution amplitude) is 
shown. As is seen the ir renormalon contributions enhance the perturbative re- 
sult and are considerable for all values of Q^, where the running coupling constant 
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method is applicable (Q^/A^ >> 1, > 2 GeV'^). The dependence of the distribu- 
tion amplitude on the scale changes only the numerical results for the ratio Rj^, 
particularly in the region > 5 GeV'^, preserving at the same time the shape of 
the curve. The ir corrections can be transferred into the scale of asiQ"^) in Eq.(15) 

Q2^^/(Q^)Q2^ f(Q^) = c^ + c,asiQ'). 

Numerical fitting allow us to get (for n/ = 3) the values of Ci, C2; 

ci ~ -8.528, C2 ~ 26.28 

The similar numerical results can be obtained also for the Rt^ using other pion 
distribution amphtudes. For the kaon a more rehable form of f{Q^) is [8] 

f{Q') = ci + C2as{Q') + C3aUQ') 
where the coefficients found after numerical fittings are 

ci ~ -1.304, C2 ~ -35.604, C3 ~ 127.25 

The results of numerical calculation for the kaon are depicted in Fig. 2. 

The obtained results have important consequences for a fate of the hard- scat- 
tering model based on pQCD at intermediate (moderate) energies (2 GeV"^ < < 
10 GeV^). It is well known that Eq.(l) was derived in the framework of pQCD 
by neglecting transverse momenta of meson's constituents and higher Fock compo- 
nents of the distribution amplitude [15]. In Ref.[17] a modified perturbative QCD 
approach was proposed, which takes into account the partonic transverse momenta, 
as well as the Sudakov corrections. In the context of this approach the pion electro- 
magnetic form factor was reexamined in works [18,19]. In these papers the authors 
used model light-cone wave functions and included also the imconventional helicity 
components hi + h2 = ±1 of the pion wave function (see, Ref.[19]). The authors 
of Ref.[18] made conclusion that the perturbative predictions are smaller than the 
experimental data, or that adding the hi + h2 — ±1 components suppresses the 
hard-scattering significantly and non-perturbative contributions will dominate in 
the present experimentally-accessible energy region [19]. But in these works in cal- 
culations the frozen coupling approximation was applied. We think that the running 
coupling constant method may improve the situation with the modified pQCD in 
the region of moderate Q^, because the ir renormalon effects work in opposite di- 
rection than factors considered in Refs. [18,19]. Therefore, we do not agree with 
Ref. [19], where the authors rule out the hard-scattering model for describing the 
light mesons form factors at moderate Q^. 

Our results do not allow us also to exclude the CZ distribution amphtude as 
the reliable model function, because in deriving of Eqs.(12),(13) we have neglected 
the second term in Eq.(lO) (Contopanagos, Sterman, Ref. [3]), as well as the second 
order corrections to Th [20]; before comparison with the experimental data [21], 
these points must be clarified. 
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APPENDIX 

The functions m„(Q^), In(Q^) for the pion have the forms (N=6), 
mi(g2) = (a + /3 + 7)% m2(g') ^ (a + 5(3 + , ui-^iQ^) = 64 (/3 + 97)' , 
m4(g2) = 16 {(3 + 147)' , m5(g2) = 23047^, mgig') = 2567^ 



and 



li(g2) = -2 (a + (3) (a + ^/?) - ^7 (58a + 90/? + 437) , 
l2(g2) = 2 (q;2 - 25/3^) - 27 (60; + 120/3 + 1357) , 
h{Q^) = 8P{a + P) + p {12a - 161/3 - 6927) , 
l4(g') = ^(3 (35/3 - a) - ^7 (14« - 417/3 + IO227) , 

k{Q') = 87 (a + 17/3 + 32I7) , l6(g') = ^7 (17177 - 3a - 20/3) . 
For the kaon we get (N=5), 

mi(g2) = {a + (3f + (7 + 5) V ^ (a + /3) (7 + 5) , 
m2(g') = (a + 5/3)' + (37 + 75)' + ^ (a + 5/3) (37 + 75) , 

QO 

m3(g') = 64/3' + 4 (7 + 95)' + y /3 (7 + 95) , 
m4(g') = 16/3' + 4005' + ^/35, m5(g') = 645', 



and 

li(g') = -2 (a + /3) (^« + ^/3) -4 (7 + 5) (7 + ^5) (9o;7 + 11q;5 + 13/37 + 15/35) , 

/ \ / 29 \ 2 

l2(g2) = 2 (a' - 25/3') - 2 (37 + 75) (7 + —5] + - {6a-f - SaS - 6O/37 - 250/35) , 

l3(g') = 8/3 (a + /3) + 2 (7 + 95) (57 - 195) + ^ (0^7 + 9«<^ + 21/37 - 67/35) , 

4 20 2 

l4(g2) = _/3 (35/3 _ + _5 (5^ + 415) + _ (-I0a5 + 10(3j + 432/35) , 
3 3 9 

l5(g') = 45 (-7 + ^5) + ^5 (a + 17/3) . 
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FIGURE CAPTION 



Fig.l The ratio i?7r for the pion (CZ amphtude) is shown as a function of . 
The sohd curves 1 and 2 correspond to -R^r with in the frozen couphng 

approximation and after the scale-setting procedure, respectively. The dashed curve 
is the ratio i?^ from Ref . [7] , where the dependence of the pion distribution amplitude 
on the hard scale was neglected. In calculations the QCD parameter A is taken 
to be equal to 0.1 GeV . 

Fig.2 The same, but for the kaon. The dashed curve is taken from Ref. [8]. 
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